Measurements of the Decay Kl e^e /i"*"^ 
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The KTeV experiment at Fermilab has isolated a total of 132 events from the rare decay Kl 
e'^e~ n'^ fj,~ , with an estimated background of 0.8 events. The branching ratio of this mode is 
determined to be (2.69±0.243tat ±0.12^;,^*) x 10"^, with a radiative cutoff of M^^^^ /M|- > 0.95. The 
first measurement using this mode of the parameter a from the D'Ambrosio, Isidori, and Portoles 
model of the Kl 7* 7* vertex yields a result of —1.59 ± 0.37, consistent with values obtained from 
other decay modes. Because of the limited statistics, no sensitivity is found to the DIP parameter 
/3. The magnitude of the angular distribution asymmetry between the e"'"e~ and i-i'^ fi~ planes, 
indicative of a CP— violating contribution to the decay, is found to be consistent with zero. We set 
a 90% C.L. upper limit of 4.12 x 10~^^ on the branching ratio of the lepton flavor-violating mode 
Kl — > e^e^/i^^^ , a factor of three improvement over the current limit from the KTeV experiment. 

PACS numbers: 13.20.Eb, 14.40.Aq, 12.15.Hh, ll.30.Er, ll.30.Fs 



The rare decay Kl ^ e+e^/i+/i^ offers the most di- 
rect means for studying the dynamics of the Kl 7* 7* 
vertex. This information is useful for models that re- 
late the Kl — > 1l~ branching ratio to n, the real 
part of the GKM matrix element Vtd H, 1 11 • This de- 
cay mode can also be used to determine the presence of 
any CP— violating contributions to the Kl 7* 7* inter- 
action Additionally, a search for the lepton flavor- 
violating counterpart Kl — > e^e^fi^fX^ provides a con- 
straint on physics beyond the Standard Model. 

In the model of D'Ambrosio, Isidori, and Portoles 
(DIP) the Kl 7* 7* form factor can be written as 

+ B ^i^i fi) 

+ ^(,?-M2)(g|-Af2)- ^ ^ 
Here, qi and (72 sue the momenta of the two virtual pho- 



tons, and Mp is the mass of the p vector meson. In this 
model, a and (3 are two arbitrary real parameters and are 
expected to be of order one. The determination of both 
a and (3 is possible through the decay Kl — > e~ fj,~ 
by examining the dilepton invariant masses and the in- 
tegrated decay rate. Knowledge of the Kl 7* 7* form 
factor is important for understanding the long distance 
contributions to Kl /^^M^ ^-nd extracting the value 
ofp[|. 

Two measurements have been made of the linear DIP 
parameter a to date, both by the KTeV collabora- 
tion. From the mode Kl j-i^ , the shape of 
the dimuon invariant mass distribution (Af^^) and the 
measured branching ratio have been used to determine 
a = -1.54 ± 0.10 111. A fit to the dielectron mass 
distribution {M^e) from Kl e^e^e^e~ determines 
a — —1.1 ± 0.6 0, where the larger error results from 
the smaller of the dielectron distribution. No mea- 
surements have yet been made of the quadratic DIP pa- 
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rameter j3. Because the effects of the /3 parameter are 
most significant in the region where both q\ and are 
large, the decay e+e^^+^^ represents the best 

means for determining (3. 

KTeV, a fixed target experiment located at Fermilab, 
collected rare decay data during run periods in 1997 and 
1999. Forty three Kl e^e^ ji^ fi'^ events were ob- 
served in the 1997 data, leading to a published branching 
ratio of (2.62 ± QAQstat ± Q-lTsyst) x 10"^ §. However, 
no attempt has been made to extract form factor infor- 
mation from this limited dataset. The results presented 
in this Letter are based on a reanalysis of the 1997 KTeV 
dataset, combined with the analysis of new data collected 
during the 1999 run. 

The two parallel Kl beams used by KTeV are cre- 
ated by focusing 800 GeV/c protons from the Tevatron 
onto a BeO target. A 65 m long vacuum region, starting 
94 m downstream from the target, defines the fiducial re- 
gion for kaon decays. Charged particles are detected with 
a spectrometer system consisting of four drift chambers 
and an analysis magnet. The hit position resolution in 
the chambers is approximately 100 ^m, while the overall 
momentum resolution is just over 1% in the range of in- 
terest. The transverse momentum kick from the magnet 
was lowered by 25% for the 1999 run for the purpose of 
increasing acceptance for four-track decay modes. 

Downstream of the spectrometer system are two trig- 
ger hodoscope planes, followed by an electromagnetic 
(EM) calorimeter. This 1.9 x 1.9 m^ array of 3100 pure 
Csl crystals has a resolution of under 1% in the energy 
range of interest. Photon vetos are located along the de- 
cay region to reject decays from particles that would miss 
the Csl calorimeter. 

Behind the calorimeter are a 10 cm thick lead wall and 
4 m of steel, the last 3 m of which serve as the neu- 
tral beam dump. A muon hodoscope (MU2) consisting 
of 56 overlapping scintillator paddles is located behind 
the dump. Following MU2 is another 1 m thick steel 
filter, behind which are two scintillator planes, one ori- 
ented horizontally and the other vertically. Known as 
MU3Y and MU3X, these planes are used for muon iden- 
tification and have 15 cm segmentation. The momentum 
threshold for muons to reach the MU3 bank has been 
measured to be 7 GeV/c. The lead and steel filters add 
up to a total of 31 hadronic interaction lengths. A more 
detailed description of the KTeV detector can be found 
elsewhere |l^ . 

The trigger requires hits in the upstream drift cham- 
bers and the trigger hodoscope planes consistent with at 
least two charged tracks. During the 1997 run, at least 
two hits were required in both MU3X and MU3Y. This 
condition was loosened during the 1999 run to allow for 
one missing hit. This change accepts events in which 
the muons are well separated in one view but happen 
to strike the same paddle in the other view. To counter 
the increased trigger rate from this change, the minimum 



number of calorimeter clusters with at least 1 GeV of en- 
ergy was raised from one in 1997 to two in 1999 ||ll|. If 
at least 0.5 GeV of energy is found in any of the photon 
vetos, the event is discarded. Events in which at least 
three tracks form a loosely defined vertex are tagged as 
candidate signal events. 

During the analysis stage, the ratio E/P is used for 
particle identification, where E is the energy deposited by 
the track in the EM calorimeter, and P is the track mo- 
mentum as measured by the spectrometer. Tracks with 
0.95 < E/P < 1.05 are identified as electrons. Tracks 
are identified as muons if they have E/P < 0.8, deposit 
less than 1.5 GeV in the calorimeter, have momentum 
greater than 7 GeV/c, and hit at least 2 out of the 3 
muon identification planes (MU2, MU3X, and MU3Y). 
Events are accepted if they contain exactly four tracks, 
with oppositely charged electron and muon pairs. 

Additional requirements imposed on the dataset in- 
clude cuts on the reconstructed kaon momentum (20 
GeV/c < Pk < 220 GeV/c) and the z position of the 
reconstructed vertex (90 m < Zytx < 158 m). Because 
the detector acceptance for decays falls off quickly 
outside of these ranges, any events observed outside of 
the boundaries are most likely misreconstructed Kl, Ks, 
or hyperon decays. To further reduce the number of mis- 
reconstructed and background events, a cut is made at 
P^ < 250 MeV^/c^, where Pt is the transverse compo- 
nent of the reconstructed kaon momentum relative to the 
kaon line of fiight. The signal mass region is defined to 
be 482 MeV/c2 < Mee^^. < 512 MeV/c^. 

Three sources of background are considered. The de- 
cay Kl ^ 7r"'"7r~7r^ (where tt^ signifies the Dalitz decay 
TT^ e+e^7) could appear as signal if the charged pions 
decay to muons in fiight, or punch through the muon filter 
and fire MU3. These events are removed from the dataset 
by cutting on an extra EM cluster in the calorimeter, in- 
dicative of the photon from the Dalitz decay of the tt". 
Two simultaneous Kl n^fi^i/ (if ^3) decays could also 
simulate signal if both pions hadronically interact and 
shower in the calorimeter, mimicking electrons. As two 
separate two-track decays are unlikely to form a good 
four-track vertex, a cut on four-track vertex quality is 
successful in removing this background. 

The most significant source of background comes from 
Kl —>■ fi'^iJ.^j events in which the photon converts to 
an pair in the material upstream of the first drift 

chamber. Monte Carlo simulations predict almost 50 of 
these events in the signal mass region at this stage of 
cuts. Requiring that the e"*" e~ hit separation at the first 
drift chamber is greater than 2 mm or that Mge > 2.75 
MeV/c^ eliminates 99% of these conversion events, while 
retaining approximately 85% of the signal (Fig. m. 

After all cuts, 132 signal events remain (Fig. 0) with 
an estimated background of 0.8 events, dominated by 
Kl — > /i~''/i~7 conversions. The background estimate is 
determined from Monte Carlo simulations. We extract 
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FIG. 1: Electron hit separation at the first drift chamber 
(in millimeters) vs. Mee (in GeV/c^) for a) Kl fi'^fi^'y 
conversion Monte Carlo , b) signal Monte Carlo, and c) data. 
Plots are logarithmic in z. The box shows the location of the 
cut used to eliminate this background. 
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FIG. 2: a) Meem^ for all data (dots) and Kl 'k'^tv~tv% 
decay /punchthrough Monte Carlo (histogram) after all cuts. 
132 events remain in the signal mass region. The vertical lines 
indicate the expected signal region, b) Close-up of the signal 
mass region. 



the branching ratio by normalizing to Kl — > tt"*" 
events, which are collected in a similar trigger. In addi- 
tion to the cuts described, for the normahzation mode the 
charged pions are required to strike the calorimeter away 
from the edges of the array (7 cm from the outer perime- 
ter, 5 cm from either beam hole), the energy deposited 
in the calorimeter by the Dalitz photon is required to be 
greater than 2 GeV, and the mass of the reconstructed 
tt" is required to be between 125 and 145 MeV/c^. Using 
these normalization events, and an acceptance of approx- 
imately 1.7% determined from Monte Carlo, the total 
yield of kaon decays within the detector from the 1997 
and 1999 runs is calculated to be (6.39 ±0.04) x 10", the 
error being purely statistical. 

Because muons are present in the signal mode but not 
in the normalization mode, accurate simulation of the ef- 
ficiency and threshold of the muon system, and a good 
understanding of multiple scattering through the muon 
filters, are crucial. These effects are studied using a com- 
bination of GEANT simulations and calibration muons 
collected with special magnet and absorber configura- 
tions. The single muon detection efficiency is measured 
to be over 99%, determined to within 0.5% of itself p^ . 

Systematic errors in the determination of the number 
of kaon decays and the signal acceptance are dominated 
by the 3.1% uncertainty in the branching ratio of the 
normalization mode. Other significant effects include a 
disagreement between data and Monte Carlo in the de- 
cay position distribution of normalization events from 
the 1999 run (1.9%), the rate of accidental activity in 
the detector (1.7%), and hmited Monte Carlo statistics 
(1.1%). Other effects include sensitivity to the vertex 
quality (1.1%), the trigger (1.0%), and the fitting proce- 
dure (0.8%). The remaining effects result from uncertain- 
ties in the background and the calibration (0.7%). The 
total systematic error on the branching ratio is 4.6%. 

Because the branching ratio and a are linked through 
the form factor, the following method is used to deter- 
mine both of these parameters. The measured branching 
ratio is determined as a function of a, where the a de- 
pendence results from the Monte Carlo acceptance cor- 
rection. Theoretical predictions for the branching ratio 
are also calculated as a function of a. Full single-loop 
QED radiative corrections are included ||l^ with a cut- 
off at Mge^^ > 0.95, and the quadratic DIP pa- 
rameter P is assumed to be zero. The intersection of 
these two curves, shown in Fig. |^, determines B{Kl — > 
e+e-ti+^i- ) = (2.69 ± 0.24,t,t ± 0.12,j,,t) x 10"^ and 
a = —1.51 ± 0.34s(at ± 0.l7syst- Varying the value of 
(3 over a reasonable range has a negligible effect on the 
results. 

Because a and (3 are connected to and g| (Eqn. |^), 
independent measurements of the form factor parameters 
can be made by studying the shape of the M^^ and Mee 
distributions of the 132 signal events. With /? fixed at 
0, the signal Monte Carlo is reweighted over a range of 
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FIG. 3: Simultaneous determination of the B(Kl — > 
e^e~ fi^ fi~ ) and the DIP a parameter. The triangles show 
values for the branching ratio, based on 132 signal events with 
0.8 background events, using different values of a to calculate 
the signal acceptance. The error bars include statistical and 
systematic errors. The dots show the theoretical dependence 
of the branching ratio on a. The intersection of the two lines 
provides the results B{Kl e+e"/.t+/i" ) = (2.69±0.24stat ± 
O.Usyst) X 10-3 and a = -1.51 ± 0.34,t„t ± 0.17,j,,t. In per- 
forming this fit, it is assumed that the DIP (3 parameter is 
0. 
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FIG. 4: a) M^^ and b) M^e (right) distributions for the 132 
signal events, compared to Monte Carlo generated with a = 
—4.53. /3 is zero for both overlays. 



values for a. Comparison to the data for each value of a 
leads to a log-likelihood distribution that is maximized 
when a = -4.53j:2:?o- A comparison between data and 
Monte Carlo at this value of a is shown in Fig. ^. Due to 
the limited statistics, we find that this analysis is insen- 
sitive to the value of the quadratic form factor parameter 

A weighted average of the two Kl — > e^e^/i+/i^ mea- 
surements of a leads to the final result a = —1.59 ± 
0.37, consistent with the two previously published val- 
ues. Combining the results from —>■ fj.'^fj.~j^ — > 
e^e~ e^e~ , and — > e^e~/i+/i^ gives a world average 
of a = —1.53 ± 0.10, a result that is dominated by the 
measurement from — > fj,'^ fi~^ . 

The Bergstrom, Masso, and Singer (BMS) model for 
the Klj*j vertex can be generalized to the two-virtual 
photon case The BMS form factor contains only 



FIG. 5: Angular distribution of the decay products for 
Kl — > e^e~ ^'^ ii~ events (dots) and Monte Carlo with no 
CP— violation (histogram) . 



one unknown parameter, ax* , which can be algebraically 
related to the DIP parameter a |5|. Using this relation, 
along with the measured value of a, we find that oik* — 
—0.19 ± 0.11 from Kl e'^ e~ fi'^ . This is consistent 
with other KTeV measurements a s well as those 

from the NA48 experiment at CERN^^. 

Some models of the Kl 7* 7* vertex allow for 
CP— violating contributions to the interaction, the pres- 
ence of which would lead to an asymmetry in the angular 
distribution of the decay products Q. The distribution 
of events in sin cj) cos (j> is shown in Fig. ^. cj) is the an- 
gle between the normals to the electron and muon de- 
cay planes in the kaon rest frame. The asymmetry A 
is defined by the ratio (7V+ - N_)/{N+ + N_), where 
Af+ (A^-) is the acceptance corrected number of signal 
events in the positive (negative) region of Fig. ||. It is 
found that A ~ — 5.3± 12.3%, which translates to a 90% 
confidence limit of |^[ < 25.5%. We conclude that no 
evidence currently exists for a CP— violating component 
of the Kl 7* 7* interaction. 

By changing the charge requirement on the lepton 
pairs, a search can be performed for the lepton flavor- 
violating decay Kl — > e^e^/i^^^ . After imposing a 
set of analysis cuts otherwise identical to those described 
earlier, no events remain in the signal region. Four-body 
phase space Monte Carlo was generated to calculate an 
overall acceptance for this mode of approximately 9%. 
This leads to a 90% C.L. limit oiB{KL e±e±/i=F/i=F ) < 
4.12 X 10"^"'^, a factor of three improvement over the pre- 
viously published limit 

In conclusion, 132 Kl e+e^/i+/i^ events have 
been observed from the 1997 and 1999 runs of the 
KTeV experiment, with an estimated background of 0.8 
events. This signal leads to a measured branching ratio 
of (2.69±0.24stat±0.12syst) x 10"^, consistent with pre- 
viously published results. In the first measurement of the 
DIP form factor parameter a using this decay mode, we 
find that a = —1.59 ± 0.37, in agreement with measure- 
ments from other modes. Because of the limited statis- 
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tices, no sensitivity is found to the quadratic DIP param- 
eter (3. In order to reduce the error on /3 to be on the 
order of one, we estimate that approximately 1,000 times 
more data would be required. Furthermore, no evidence 
is found for CP— violating contributions to the Kl 7* 7* 
interaction. Finally, we constrain the branching ratio of 
the Icpton flavor-violating mode Kl e^e^/i^/i^ to 
< 4.12 X IQ-", at 90% confidence. 
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